Mon. Not. R. Astron. Soc. 000, 000-000 (2001) Printed 1 February 2008 (MN KTpX style file vl.4) 



U 9 B and r band luminosity functions of galaxies in the Coma 
cluster* 



o 
o 

a* 

oo 

(N 

(N 
> 

m 
o 



a: 

6 

CO 



Marco Beijersbergen/t Henk Hoekstra, 2 3 Pieter G. van Dokkum, 4 i 
and Thijs (J.M.) van der Hulst 1 

1 Kapteyn Astronomical Institute, P. O. Box 800, 9700 AV Groningen, The Netherlands 
2 CITA, 60 St. George Street, Toronto, M5S 3H8, Canada 

3 'Department of Astronomy, University of Toronto, 60 St. George Street, Toronto, M5S 3H8, Canada 

4 California Institute of Technology, Mail Stop 105-24, 1200 E California Blvd Pasadena, CA 91125, USA 



Version 1 February 2008 



ABSTRACT 

We present a deep multi-colour CCD mosaic of the Coma cluster (Abell 1656), covering 5.2 
deg 2 in the B and r bands, and 1.3 deg 2 in the U band. This large, homogeneous data set 
provides a valuable low redshift comparison sample for studies of galaxies in distant clusters. 
In this paper we present our survey, and study the dependence of the galaxy luminosity func- 
tion (LF) on passband and radial distance from the cluster centre. The U, B and r band LFs 
of the complete sample cannot be represented by single Schechter functions. For the central 
l ioo kP c ' we fi 11 ^ best-fitting Schechter parameters of M^j 

-,+0.018 

-0.028 



area, r < 245 
a v = -1.32^ 



M* = 



-19.79±g;if and a B = 



mmtoll and 



1.37±°;°* and M* = 



and a r = — 1.161o'q^ 9 . The LF becomes steeper at larger radial distance from the cluster cen- 
tre. The effect is most pronounced in the U band. This result is consistent with the presence 
of a star forming dwarf population at large distance from the cluster centre, which may be in 
the process of being accreted by the cluster. The shapes of the LFs of the NGC 4839 group 
support a scenario in which the group has already passed through the centre. 
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1 INTRODUCTION 

Clusters of galaxies are important laboratories for studies of galaxy 
evolution. The galaxy population in clusters is very different from 
the population in the field, suggesting that galaxy formation and 
evolution are a strong function of the environment. 

The galaxy luminosity function (LF) should be an excellent 
tracer of environmental effects. Knowledge of the shape of the LF 
(in different bands) is a powerful tool for studies of galaxy evolu- 
tion. Specifically, one can look for correlations between the shape 
of the general LF and environmental or cluster properties. In gen- 
eral the LF drops steeply at bright magnitudes and rises gradually at 
faint er magnitudes as described by a Schechter function (S chechter 
1976). Sometimes, however, features such as bumps and a steeply 
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rising faint part are found, which cannot be adequately fitted by a 
single Schechter function. Galaxies in the cluster core region are 
expected to have a different merger history than the galaxies popu- 
lating the cluster outskirts where it blends into the field. This should 
be reflected in differently shaped LFs for dense and less dense re- 
gions within a cluster. Lopez-Cruz et al. (1997) propose that the 
flat faint end slopes found in rich c luster s result from the disruption 
of dwarf galaxies. Biviano et al. (1995) report a dip in the bright 



part of the general LF for rich cluster s whic h is not seen in LFs of 
poor clusters or in the field. Andreon (1998) verified the invariance 



of the shape of the bright part of the type-dependent LF in a large 
range of environments from the field to the cores of clusters several 
orders of magnitude denser. The determination of the exact shapes 
of LFs is difficult as the faint ends suffer from uncertainties in the 
contamination by field galaxies. 

The Coma cluster (z = 0.023, richness class 2) is the richest 
of the nearby clusters and ideal to study environmental effects (for 
a detailed overview of research on the Coma cluster see Mazure 
et al. 1998). Previous large field studies of the photometric proper- 
ties of the galaxies in the Coma cluster have been based on pho- 
tographic plates (e.g. Godwin, Metcalfe & Peach 1983; Lugger 
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Figure 1. Layout of the four science CCDs and the autoguider CCD as used 
in our survey of the Coma cluster (rotator angle at 270°). The four science 
CCDs cover an area of 34' X 34'. 



1989). CCD studies yield better photometric precision, but have 
hitherto been limited to relatively small areas, mainly focused on 
the central regions (e.g. Thompson & Gregory 1993; Biviano et al. 
1995; Bernstein et al. 1995; Lopez-Cruz et al. 1997; Seeker et al. 
1997; Lobo et al. 1997; de Propris et al. 1998; Trentham 1998a, 
1998b; Mobasher & Trentham 1998; Andreon 1999). We combine 
the benefits of wide field coverage with the photometric accuracy 
attainable with CCDs, by using the Wide Field Camera (WFC) on 
the Isaac Newton Telescope (INT) to image a large part of the Coma 
cluster from the dense core to the outskirts of the cluster. Our data 
offers the interesting possibility to study the LF for a large range of 
environments within the cluster. The primary goals of our survey 
are to study the LF, the colour-magnitude relation and to comple- 
ment a large HI survey conducted with the Westerbork Synthesis 
Radio Telescope (WSRT). 

In this paper we present our data and use our catalogue to con- 
struct U, B and r band LFs for various regions of the Coma clus- 
ter area. Our study complements previous studies of the LF of the 
Coma cluster by providing accurate, deep CCD photometry in three 
bands, covering an area of 5.2 deg 2 . Our wide field U band data is 
of particular importance, because the U — R colour is a sensitive 
indicator of the presence of young stellar populations. Furthermore, 
the rest frame U band remains in the optical window out to z ~ 1.3, 
which makes the U band observations very useful for comparison 
to high redshift clusters. Previous U band studies of Coma were 
limited to small samples of individual galaxies ( Bower, Lucey & 
Ellis 1992). Our mosaic has similar resolution (in kpc) and size (in 
Mpc) as large HST WFPC2 mosaics of high redshift clusters (e.g. 
van Dokkum et al. 1998) and provides a valuable low z zeropoint 
for studies of the evolution of galaxy morphology, colour and lu- 
minosity. For the first time a U band LF for such a large area of the 
Coma cluster is determined. In addition we study the dependence of 
the LF on passband and projected distance from the cluster centre. 

Throughout this paper we use Ho = 100 ftioo km s _1 Mpc -1 
and coordinates are given in equinox J2000.0. 



2 OBSERVATIONS 

The data were collected during the nights of March 19-22 1999 
with the Wide Field Camera (WFC) on the Isaac Newton Telescope 
(INT), on Rogue de Los Muchachos on the island of La Palma 
(Spain). The WFC consists of four 2048 x 4100 pixels EEV CCDs 
and a fifth CCD for autoguiding. The layout of the chips is shown 
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Figure 2. Layout of the 25 overlapping pointings covering an area of ~ 5.2 
deg 2 . Thick lines delineate pointings which were observed in the U band, 
as well as B and r. Pointing 11 served as reference field for the photo- 
metric calibration. The position of NGC 4839 is indicated by the square. 
Coordinates are given relative to a = 12 h 59 m 43 s , S = +27°58'14". 



in Fig.Q. It is designed to provide a large field survey capability for 
the prime focus of the INT. The sky coverage is 4 x 259' 2 with a 
plate scale of 0'.'333 pixeP 1 . The camera covers 80 per cent of the 
unvignetted field of the INT. 

We have imaged a mosaic of 25 overlapping pointings cover- 
ing a total area of ~ 5.2 deg 2 or 2.8 x 2.8 h j qq Mpc. Broadband 
filters RGO U, Harris B and Sloan r were used for the 6 point- 
ings covering the core and the south-west group (containing NGC 
4839). The other 19 pointings were observed with Harris B and 
Sloan r. For each pointing two exposures, offset by ~ 1', were 
taken. This ensured that we would be able to determine the rela- 
tive zeropoint offsets between exposures and chips and there would 
be very few gaps in the final mosaic. Furthermore, galaxies in the 
overlapping regions give us good estimates of the errors due to pho- 
ton noise and flatfielding. Integration times were 2 x 300 s in B, 
2 x 600 s in r and 2 x 900 s in U. The layout of the field is shown 
in Fig. ^| Thick lines delineate pointings which were observed in 
the U band, as well as B and r. 

The overall weather was good with an average seeing during 
the first two nights of ~ 1.4" and ~ 1.9" during the last two nights 
in the r band. The strategy was to observe the central parts of the 
cluster in U, B and r under photo metric conditions. During the first 
night the standard Landolt (Il992j) fields Sal07-602 and Sal01-427 



were also observed at similar airmasses as the central fields. These 
standard fields contain many standard stars and give a good spread 
over the four CCD chips. This is important for the absolute calibra- 
tion, since no two chips have exactly the same characteristics. In 
fact it is known that the WFC is non-linear. 

The Coma fields were observed when the cluster had risen to 
an airmass < 1.8. At the start of nights 3 and 4 we took exposures 
of an empty field (see Section 4.4). This field was later used to 
correct U band galaxy counts for foreground/background contami- 
nation. 
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3 INITIAL REDUCTION 

3.1 Flatflelding 

Prior to any reduction the images were inspected to judge their 
quality. For each chip the bias frames per night were checked for 
repeating structures and count levels. The bias-structures were also 
visible in the science images. We combined all bias frames to cre- 
ate an average master bias frame for each night and each chip. Af- 
ter several tests we decided that the bias subtraction had to be per- 
formed in two steps; we first subtracted the overscan level by fitting 
a second order polynomial to the columns and then subtracted the 
(overscan subtracted) master bias frames. The second step success- 
fully removed the remaining bias structures for chips 1,2 and 4. 
Chip 3 suffered from a few bad columns. These were removed by 
interpolation or by aligning and combining the images taken with 
a small offset. The rms noise level in the master bias frame was 
negligible compared to the rms noise in the science frames. 

Accurate flat field images were created from the science 
frames. We have obtained 52 science images of the Coma area per 
chip per band for the B and r bands. For the U band we have 23 sci- 
ence images, including the empty fields. We constructed a B and r 
band masterflat by scaling these images by the mode, rejecting the 
lowest 10 and highest 15 pixelvalues, followed by median filtering. 
The U band masterfiats were constructed in the same way, but there 
we rejected the lowest 5 and highest 10 pixelvalues before median 
filtering. 

As a test we compared our flatfields to the ones of the Wide 
Field Imaging Survey which are available from the WFC archive^ 
Differences between the flat fields were at the level of 10~ 4 or less, 
assuring that we have good quality flats. 

After flatfielding, approximately 2.5 per cent of the area of 
chip 3 still suffers from ~ 3 per cent variations in the background, 
because this chip is partly vignetted. Because of the overlap be- 
tween neighbouring pointings this has very little impact on our fi- 
nal mosaic. On average, the images are very flat with less than 0.8 
per cent variations in the background. 

3.2 Photometric calibration 

The ultimate aim was to obtain a homogeneous photometric scale 
over the whole area covered. We dedicated the photometric anchor 
point to be chip 4 of pointing 11 (Fig. ^) at airmass 1.000. This 
field was observed during the first night under photometric condi- 
tions at a time close to when the standards were taken. First, we 
determined the photometric offsets between the four chips of the 
WFC, using objects in overlapping regions between the two offset 
exposures for each pointing. For the U band the number of bright 
objects in overlapping areas was low. Therefore, we could not use 
objects to determine the offsets, since this would have given unreli- 
able results. Instead, we used the sky level to determine photometric 
offsets. First, we made sky-images by filtering out all objects of the 
science images in the same way as for the flatfields. Comparison of 
the mean skylevels of the sky images gave the offsets. The result- 
ing offsets are listed in Table [lj Errors indicate the uncertainty in 
the mean of all independent measurements. The la spread between 
independent measurements of pairs of different pointings is much 
larger, and is caused by varying conditions between observations 
of neighbouring pointings. The offsets depend on the passband, be- 
cause the response curves of the four CCDs are not identical. 

§ http://archive.ast.cam.ac.uk/wfsurvey/wfsurvey.html 



Table 1. Photometric offsets relative to chip 4 

chip U B±cr r±<r 



+0.347±0.015 
+0.170±0.015 
+0.O16±0.015 



+0.337±0.013 
+0.410±0.013 
+0.325±0.014 



+0.446±0.011 
+0.501±0.017 
+0.376±0.011 



Table 2. Average extinction coefficients for La Palma 

band A c g (A) extinction/airmass (mag/airmass) 



u 


3610 


0.45 


B 


4361 


0.20 


r 


6216 


0.08 



We applied these offsets to bring all chips to the same scale. 
Only the first two nights were photometric and we proceeded by 
determining the photometric offsets between neighbouring point- 
ings. The mean offset values are based on the 8 (or more) brightest 
objects in the overlapping areas in order to guarantee a reliable de- 
termination. We found differences of a few hundredths of a magni- 
tude between fields taken at significantly different times (airmasses, 
conditions) in B and r up to more than a tenth of a magnitude in U. 
We mapped all the photometric offsets per pointing per band and 
scaled all pointings in such a way that the errors were spread over 
the total observed region and did not accumulate towards the edges 
of the covered area. The final result is that all B and r band images 
have the same zero-point to ~ 0.04 mag and the U band images 
to ~ 0.06 mag. This is shown for the r and U bands in the right 
panels of Fig. pi 



The photometric calibration was done using Landolt (1992) 
standard stars. We did not have sufficient standard stars on all the 
chips to solve for the extinction coefficients so we set the extinction 
coefficients to constant values. We used average extinction coeffi- 
cients at the effective wavelengths of the filter+CCD system, as 
listed in Table |^, to derive the colour terms and zero-points. The 
zero-points are listed in Table [| The colour terms were found to be 
very small (~ 0.01) and are neglected. Since Coma lies close to the 
galactic pole the extinction is insignificant for our purposes (0.043, 
0.034 and 0.021 mag in the (J, B and r bands respectively). Hence, 
we did not apply extinction corrections. 



4 OBJECT CATALOGUES 

The data volume is considerable and requires a fully automated 
pipeline for data handling. After the standard reduction steps the 
processed images were presented to software implemented at the 
Leiden Data Analys is Center (L DAC) for pipeline processing of 
overlapping images (Deul 1998). The processing software is a se- 
ries of programmes that, when run in a chain, derive source param- 
eters for any given set of input frames of a given passband. Only the 
first step is an actual data reduction task in that it reduces the infor- 
mation to be processed from image data to catalogue data (an object 
list). All but the last of the following pipeline routines add informa- 



Table 3. Zero-points (ZP) of chip 4 for night 1 

ZP±CT 



U 
B 



23.157±0.026 
24.918±0.031 
24.822±0.073 
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Figure 3. Astrometric and magnitude scatter in the r and U bands for objects in overlapping regions. The circle denotes the astrometric accuracy cr r ~ 0.3' 



tion to the catalogue(s) increasing its size and information content. 
The final step merges the catalogues to create a multi-colour object 
catalogue with astronomically meaningful source information. The 
pipeline processing is performed on the full data set. Along the way 
we must decide on values for various parameters. We describe the 
pipeline steps in more detail below. 



4.1 Detection of objects 

We used SExtractor to automatically detect objects on each chip 
individually. The complete analysis of an image is done in six 
steps. First, a model of the sky background is built and parame- 
ters describing the global statistics are estimated. Then the image 
is background-subtracted, filtered and thresholded. Detections are 
then deblended, cleaned, photometered, classified and written to 
the final catalogue. For specific details o f each of these steps the 
reader is referred to Bertin & Arnouts (1990 or the SExtractor 
user's guide. 



Objects were identified as the peaks in the (background sub- 
tracted) convolved images that were higher than a given threshold 
above the local background. We used most of the standard SEx- 
tractor settings except for the memory parameters which had to be 
increased in order to detect the large cD galaxies successfully. The 
seeing parameter had to be adjusted for each image individually 
for a good star/galaxy separation. For all objects positions, magni- 
tudes, basic shape parameters and star/galaxy (S/G) classifiers were 
determined and written to a catalogue. SExtractor produced a cata- 
logue for each image which was subsequently converted to a format 
suitable for pipeline reduction. As a first pipeline processing step, 
all information of the individual image source extractions plus all 
original FITS image header information is combined into one sin- 
gle output catalogue per band. 



4.2 Astrometry 

Astrometric calibration was performed by pairing the input posi- 
tion catalogue (USNO-A2) with the extracted object information. 
For multiple band processing intercolour pairing is done first and 
between frames in overlap, overlap pairing is performed as well. 
The derived astrometric solution is then applied to all the objects 
in the set of frames. Sky coordinates (RA and Dec), as well as cor- 
rected geometric parameters are calculated. The final precision of 
this calibration depends on the accuracy of the source extractions, 



input catalogue accuracies and the correctness of the functional de- 
scription of the distortions. We find cr r ~ 0.3" as shown in the left 
panels of Fig. ^ 



4.3 Final object catalogue 

The final step in creating a catalogue containing useful multi-colour 
astronomical data is to merge all catalogues to get one source per 
position on the sky. All information referring to the same astro- 
nomical object is gathered and merged. Position information is 
a weighted mean of all detections where the weighting is based 
on detection signal-to-noise and detection environment conditions. 
Details of the merging are configurable. When the brightness con- 
trast between overlapping galaxies is too low the software is not 
able to deblend them correctly, resulting in an erroneous catalogue 
entry. On the other hand, setting the deblend contrast parameter too 
low causes SExtractor to consider bright star forming regions in a 
galaxy as separate objects. Then, there is the possibility that de- 
blended objects are merged by the pipeline software when creating 
the final object catalogue. This can happen in cases where objects 
are small compared to the errors in position and shape parameters. 
We carefully inspected dense regions in our mosaic and conclude 
that erroneous catalogue entries are rare and consequently do not 
affect our results. 

Star/galaxy separation was performed based on SExtractor's 
stellarity index. For bright objects stars and galaxies are easily sep- 
arated, but towards fainter magnitudes and for bad seeing the divi- 
sion is not so clear. For most purposes we consider the objects with 
a S/G classifier value smaller than 0.8 in the r band (best seeing) 
to be galaxies. Bright (saturated) stars tend to have a S/G classifier 
smaller than this. We filtered these out by demanding a S/G classi- 
fier < 0.1 for the brightest magnitudes. We visually inspected all 
bright objects to conclude that there is no contamination by stars up 
to at least m r = 15. Beyond m ~ 19 mag we probably still have 
a fraction of stars in the sample, but for most purposes it is better 
to be contaminated by a (small) fraction of stars than to reject com- 
pact galaxies. Histograms of the catalogue's raw number counts are 
shown in the panels of Fig. Q Troughout this paper we use SExtrac- 
tor's MAG_BEST as magnitude estimator. The final mosaic, com- 
posed of all pointings observed in the B and r bands, is shown in 

Kg. a 
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Figure 4. Histograms of the normalized raw r, B and U band number counts. N tota i gives the total number of catalogue entries for each passband. 



4.4 Control fields 

For a survey as large and deep as ours, foreground/background sub- 
traction can only be treated statistically, since only the brightest 
galaxies have measured redshifts. The main source of uncertainty 
in the determination of LFs comes from number statistics and back- 
ground variance. Usually, flanking fields are used to get an estimate 
of the foreground/background correction and its variance. We ex- 
pect that over the large area of the WFC the effects of cosmic vari- 
ance are not important. 

For the U band we have used 4 x 900 s observations of an 
empty field located at a = 8 h 00 m 00 s , S = +50°00'00" span- 
ning ~ 980 arcmin 2 . For both the B and r bands we have made use 
of the Wide Field Survey (WFS) data archive to get images of ran- 
dom fields. The B band control images are eight 600 s exposures, 
spanning ~ 2000 arcmin 2 and centered at a = 12 fe 55 m 55 s , 5 = 
+27°01'40" and a = 12 h 53 m 40 s , 8 = +26°20'00". The r band 
control images are four 600 s exposures, spanning ~ 1000 arcmin 2 
centered at a = 16 fe 04 m 26 s , 5 = +54°49'59". For the B and r 
band control fields we rely on the (photometric) reduction of the 
four chips by the WFS team. 

All control images were pipeline reduced to give control 
catalogues which were filtered with the same criteria as the 
Coma fields. We applied relative extinction corrections ( Schlegel, 
Finkbeiner & Davis 1998) to bring the extinction in the control 
fields into agreement with the Coma extinction, not zero extinction. 

4.5 Completeness 

The histograms shown in Fig. ^ suggest that the limiting magnitude 
is ~ 22.5 mag for all bands. However, due to less than perfect con- 
ditions in the last two nights this limit does not apply to the total ob- 
served area. To have uniform completeness for the total mosaic we 
compute LFs up to mu m = 21.73. The positions of these objects, 
relative to the cluster centre, are plotted in Fig. |6| The cluster is vis- 
ible as a density enhancement on top of a uniform background. We 
have only considered clean detections (with good extraction flags) 
with more than 5 connected pixels above the background threshold. 

Detection limits in general depend on the interplay between 
scale length/effective radius, magnitude and inclination, i.e. sur- 



Table 4. Total cluster counts and field counts estimates per deg 2 
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Figure 6. Distribution of objects with me < 21.73. Coordinates are g' 
relative to a = 12 h 59 m 43 s , 6 = +27°58'14". 
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Figure 5. Overview of the total area covered in the B and r band with annuli overlaid. Radii are: I: 0-0.2, II: 0.2-0.3, III: 0.3-0.42, IV: 0.42-0.74 and V: 0.74-1.4 
degrees (1° = 1.22 Mpc at Coma distance). The circular area centered on NGC 4839 has a radius of 0.15°. 



face brightness. Edge-on galaxies of a certain magnitude are easier 
to detect than their face-on counterparts. Deep surveys have shown 
the existence of galaxies with surface brightnesses fainter than the 
night sky. These low surface brightness (LSB) galaxies are more 
likely to be missed at a given magnitude and inclination than their 
high surface brightness counterparts. Most of the LSB galaxies in- 
vestigated in any detail are either late-type and disk dominated (de 
Blok et al. 1995; McGaugh & Bothun 1994), or giant, Malin-l-like 
galaxies (Sprayberry et al. 1995; Pickering et al. 1997). 

In order to estimate these effects we followed two approaches. 
First we generated a r band luminosity-surface brightness plot for 
all objects (mostly galaxies) in our catalogue. The data in our cat- 
alogue do not contain estimates of the central surface brightnesses, 
but we can compute mean surface brightnesses for all objects based 
on total flux and area above the analysis threshold. In Fig. H we 



show the resulting luminosity-surface brightness plot along with 
the relevant magnitude and surface brightness selection boundaries 
we used. The dashed lines represent the luminosity-surface bright- 
ness relations for objects with the minimum detectable size (the 
seeing disk). We show minimum size detection lines for a seeing 
of 1.5" and 1", from left to right, and it can be seen that occasion- 
ally the seeing drops even below 1". The mean isophotal limits for 
the U and B bands are 25.10 and 25.42 mag arcsec~ 2 respectively. 
From this figure it is evident that we are not likely to miss galax- 
ies due to their low surface brightness. Second we generated arti- 
ficial elhpjical_gala)de^_wath_^^ Vaucouleurs 
1948) light profiles and spiral galaxies with exponential light pro- 
files. We then added these galaxies into empty regions of our Coma 
images and verified whether SExtractor could recover these using 
the same selection criteria as for Coma fields. From our simulations 
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Figure 7. Luminosity versus mean surface brightness plot. The selection 
boundaries are indicated by the solid lines. The dashed lines represent the 
minimum size detection lines for a seeing of 1.5" and 1" (from left to 
right). 



Table 5. Best-fitting Schechter parameters for the luminosity functions 
of the complete sample 



Filter 



M* 



n/ifooMpc^mag- 1 ] x 2 



R _1QOQ+ - 36 _1 ^9+0.056 

D ry-U»_o.40 1-0 -0.049 



r _90 R'3+0-26 i 99 +0.034 

r zu.oo_ 34 i.zz_ Q 03g 



11.7 ± 1.05 27.8 
13.7 ±0.02 5.35 
20.44 ±0.48 2.87 



we estimate that we are able to detect ellipticals with effective radii 
r c ~ 3 h^Q kpc or r c ~ 9" down to ~ 19.5, ~ 19.5 and ~ 19 mag 
in the U, B and r band respectively. Furthermore, we are able to de- 
tect dwarf galaxies modelled as exponential disks with hd ~ 1 kpc 
down to ~ 19.5, ~ 20 and ~ 20 mag or ^ = 23.8, 24.3 and 24.3 
mag arcsec -2 in the U, B and r band respectively. At the distance 
of Coma, m — 20 corresponds to M = —14.2. Disk dominated 
LSB galaxies typically have /id ~ 3 kpc and /xo,b ~ 23. 2 mag 
arcsec -2 . Assuming typical colours as in the de Blok et al. (1995) 
sample, this corresponds to /Lto,u ~ 23.1 and fio tt ~ 22.4 mag 
arcsec -2 . Our simulations show that such LSB galaxies are within 
our detection limits. 

To illustrate the data quality we show examples of galaxies in 
Fig. ^[ These have been drawn from our catalogue for magnitude 
bins Mb= -21 to -13 in steps of one magnitude, two examples per 
bin. 



5 TOTAL LUMINOSITY FUNCTIONS 

We constructed foreground/background corrected LFs by subtract- 
ing counts in the control catalogues from the Coma counts. From 
Fig. ^ it is obvious that the number of foreground/background 



Table 6. Schechter parameters for the central luminosity functions 

Filter M* a (f>*[h^ 00 Mpc -2 mag -1 ] x 2 



U 
B 



-18.60 



+0.13 



+0.018 
-0.028 

-0.024 
-0.016 

90 C7+0-12 i 1K +0.012 



-19.79 



+0.18 



-1.32 



-1.37" 



135.1 ±4.3 1.78 
61.7 ±2.4 1.71 
130.3 ±11.6 1.34 



counts is a strongly varying function of magnitude. Especially the 
faintest points of the LFs are largely influenced by an inaccurate de- 
termination of the amount of contaminating galaxies. We used very 
large control fields, and large bin-widths compared to the photo- 
metric uncertainties. Therefore, effects caused by errors in the de- 
termination of foreground/background counts will be suppressed. 
Furthermore, we carefully calculated effective areas for all control 
fields and the Coma mosaic by counting unblotted pixels on each 
frame. Pixels which reduced the area for object detection, e.g. dead 
columns, saturated stars etc., were blotted prior to the pixelcount- 
ing. The total cluster counts and the field counts estimates are listed 
in Table fij 

In Fig. ^ the LFs for the complete data set are shown. Solid 
lines represent best-fitting Schechter functions with parameters as 
listed in Table ^. The faint end slopes increase towards shorter 
wavelengths. The best-fitting Schechter functions are rather poor 
representations of the data for all bands, with several points lying 
more than la away from the best-fitting values. The x 2 statistic 
confirms the eye's impression that the LFs for the complete data 
set cannot be represented by single Schechter functions. These are 
the first accurate determinations of LFs for such large areas of a 
cluster. Furthermore, to our knowledge the U band LFs presented 
here are the first ever published for the Coma cluster. 



6 DEPENDENCE OF LUMINOSITY FUNCTIONS ON 
RADIAL DISTANCE FROM THE CLUSTER CENTRE 

In order to study the dependence of the LF on radial distance from 
the cluster centre we defined five areas as annuli with varying 
widths and radii projected on the cluster centre, as shown in Fig.pl 
For each of these areas we carefully determined the effective area 
for galaxy detection and constructed the corresponding LFs. The 
contamination by foreground/background galaxies becomes severe 
towards the outskirts of the cluster. Their numbers become equal or 
greater than the Coma counts at r > 19, r > 18 and r > 17 in 
annulus I, III and IV, respectively. 

We arbitrarily defined the core of the Coma cluster as the area 
with r < 245 ft.]~ 00 kpc (area I). This is comparable in size to 
the total observed areas in previous studies. In Fig. [uj we show 
the central LFs with best-fitting Schechter functions overplotted as 
solid lines. The insets show the 68, 95 and 99 per cent confidence 
levels for the M* and a parameters resulting from the fit to the 
binned data. The best-fitting parameters are listed in Table ^. The 
bright end of the B band central LF is not adequately represented 
by the Schechter function. The faint end, however, is well repre- 
sented by the Schechter fit. We do not c onfirm the dip, located at 
Mb ~ —17.2, reported by Biviano et al. ( 1995 ). We stress that our 
LF has been determined with completely different data and meth- 
ods, complicating any direct comparison of results. Their estimate 
of the faint end slope («b = — 1.3±0.1) is, however, in agreement 
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Figure 9. LFs for the total area observed. The solid lines represent best-fitting Schechter functions with parameters as given in Table ^[ The dashed lines 
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Figure 10. LFs for area I. Solid lines represent the best-fitting Schechter functions with parameters as listed in Table ^[ The insets show the 1, 2 and 3cr contour 
levels of the best fitting Schechter function parameters. N tota i gives the estimated number of Coma galaxies up to -15.2. 



with our result. Comparison of our B band LF with t he typical rich- 
ness class 2 composite cluster LF ( Trentham 1998b| ) shows that the 
faint ends are consistent up to the completeness limit. Beyond the 
completeness limit the composite LF rises more steeply than the LF 
we have derived. Our value of the fa int en d slope of the central r 
band LF is in agreement with Lugger ( 1989) (qr ~ — 1.19±0.17), 
but somewhat s hallow er than the slopes d erived for Coma by e.g. 
Bernst ein et al. (1995), Lopez-Cruz et al. (1997) and Seeker et al. 
(1997) who all find lvr ~ —1.4. We stress that a direct compari- 



son is hampered by the fact that these authors use different areas or 



composite LFs. To demonstrate more clearly that th e LF is mainly 
shaped by the area we extracted Trentham's ( 1998a ) area and con- 
structed B and r band LFs. In Fig. [ll] we compare the results of 
both determinations of the LFs. Up to the completeness limits the 
LFs are very similar showing that not the methodology employed, 
but the area chosen has the most crucial effect on the shape of the 
LF. 

In general a single Schechter function is a reasonable repre- 
sentation of all the central LFs with 2 points lying lcr or more from 
the best-fitting value for all bands. It is expected from LF studies of 
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Figure 11. Comparison of our data to those of Trentham ( 1998a ), restricting our data to the same area covered by the Trentham ( 1998a ) survey. There is very 



good agreement between the two independent determinations of the Coma LF in this small region. 



composite dense and loose clusters and of single rich and poor clus- 
ters that the LF shape depends on environment. Below, we study the 
change in the shape of the LF as function of position in the cluster. 
We will examine whether differences can be attributed to effects of 
the local environment. 

In the panels of Fig. [l2| the LFs corresponding to the annuli II 
to V of Fig. ^| are shown for all bands. It is clear that these LFs are 
not simply scaled versions of the central LFs. 

The U band LFs are sensitive to star forming galaxies, and 
are therefore a poor indicator of the underlying mass distribution. 
They show significant curvature and at some radii dips, as reported 
for other bands (e.g. Biviano et al. 1995; Andreon & Pello 1999). 
Because of these dips these LFs cannot be adequately fitted by sin- 
gle Schechter functions. 

The B band LF of annulus II still resembles the central LF, 
but in annuli III and IV the faint end behaves differently. Even fur- 
ther out, in annulus V, the galaxies with Mb ~ —19 or brighter 
become very rare and the faint end seems to steepen again, but only 
marginally. 

The r band LF of annulus II is relatively flat. At larger radial 
distances from the cluster centre the LFs become much steeper. 

We quantified these trends as follows. We fitted a power law 
function (b 10 aM ) to the faint ends of the LFs in order to study 
their dependence on radial distance from the cluster centre. The LFs 
were fitted for Mu > -18, Mb > -19 and M T > -20, respectively. 
In Fig. [l3| we plot the power law slopes as function of cluster radius. 
We omitted the slope of the B band LF for area IV: the value of the 
slope is extremely sensitive to the fitting region, and hence not well 
constrained. In general, the faint end slopes become steeper towards 
larger cluster radii. 




0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 
Radius [h^ ' Mpc] 

Figure 13. Power law slopes a as function of cluster radius. The faint ends 
of the LFs become somewhat steeper towards the cluster outskirts. 



7 COMPARISON WITH FIELD LUMINOSITY 
FUNCTIONS 

In the field the galaxy density is orders of magnitudes lower than in 
the cores of clusters like Coma. Galaxies in a dense cluster environ- 
ment are likely to follow different evolutionary paths than galax- 
ies in low density fields. It is therefore interesting to investigate 
whether this is reflected in their LF shapes. LFs of the field pop- 
ulation have recently been measured by a number of surveys (e.g. 
Loveday et al. 1992; Lin et al. 1996; Lin et al. 1997; Geller et al. 
1997; Marzke et al. 1998; Lin et al. 1999). Despite the large sam- 
ples that were used to measure the LFs, controversy on the shape re- 
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Figure 12. LFs for the annuli of Fig. ^| In the top of each figure we indicate: filter, annulus number and estimated number of Coma galaxies up to -15.2. The 
lines represent the fits to the faint ends. 



mains. One of the largest local redshift samples of galaxies selected 
from CCD images is provided by the Sloan Digital Sky Survey 
(SDSS) (Blanton et al. 2001). The SDSS survey has multicolour 
v! ,g' ,r' ,i' , z' photometry, whereas most of the older surveys mea- 
sure LFs only for the B and R bands. We transformed the SDSS 
LFs to our pho tomet ric system using the transformations given in 
Fukugita et al. ( 199C ). Similarly, values for the slopes in U, B, and 



r were obtained by interpolating between the nearest SDSS filters. 
The resulting SDSS LFs are shown in Fig. ^| as dashed lines. For 
clarity the LFs were renormalized by adding 2.8 to log (/>(M) in 
all bands. The faint end slopes of the field LFs increase towards 
shorter wavelengths, although less significantly than in Coma. The 
faint end slopes of the r and B band LFs of field and Coma galax- 
ies are very similar. Interestingly, the faint end of the Coma U band 
LF appears to be steeper than that of the field LF. This may indicate 



that dwarf galaxies at Mb r 
dwarf galaxies in the field. 



-13 in Coma have similar colours as 



8 GALAXY DISTRIBUTION 

The best indicator of the underlying mass distribution in our data 
set is the r band, since it is least influenced by episodic star forma- 
tion. We use this band to investigate the p rojected galaxy density 
distributions. Following Driver et al. (1998 ) we separate our sample 



into giant and dwarf galaxies. We define giant galaxies as objects 
with —23.5 < M r < —19.5 and dwarf galaxies as objects with 
— 19.5 < Mi < —16.5. In Fig. [l4| we plot the projected distribu- 
tions of these two types on the sky. In panels (a) and (b) of Fig. [l5| 
we plot their foreground/background corrected projected densities 
and in panel (c) their ratio as function of distance from the cluster 
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Figure 14. (a) Filled hexagons represent galaxies with —23.5 < M r < 
-19.5 (giants). The 3 large hexagons represent NGC 4889, NGC 4874 and 
NGC 4839. (b) Crosses represent galaxies with -19.5 < M r < -16.5 
(dwarfs). Overplotted are annuli with radii: 0-0.2, 0.2-0.3, 0.3-0.42, 0.42- 
0.53, 0.53-0.74, 0.74-1.1 and 1.1-1.4 degrees (1° = 1.22 Mpc at 

Coma distance). 



centre. In panels (a) and (b) we also show isothermal profiles for 
comparison. The giant galaxies have a relatively small, sharply de- 
fined, area of high overdensity limited to the central cluster area. 
At a distance of 0.37 h^[ 00 Mpc from the cluster centre their density 
abruptly drops by ~ 70 per cent and then decreases continuously. 
At distances from the cluster centre larger than 0.37 h^ 00 Mpc the 
dwarf-to-giant ratio (D/G) could become somewhat larger than in 
the core of the cluster, but this is not significant given the large 
uncertainties. 
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Figure 15. (a) The number density of dwarf galaxies as function of dis- 
tance from the cluster centre, (b) The number density of giant galaxies as 
function of distance from the cluster centre, (c) The dwarf-to-giant ratio 
(D/G) as function of distance from the cluster centre. The NGC 4839 group 
is indicated by a square. Dashed lines correspond to isothermal profiles. 
Horizontal errorbars indicate bin widths. 



8.1 NGC 4839 group 

It has long been known that at ~ 40' south-west of the cluster cen- 
tre a secondary concentration of galaxies exists (e.g. Wolf 1901). 
This is the group associated with the cD galaxy NGC 4839. The 
presence of this group is clearly visible in Fig. fla: the group's den- 
sity lies well above the average density at that cluster radius and 
is comparable to the central densities. We used a circular area of 
255 arcmin 2 centered on NGC 4839 (Fig. to determine LFs. The 
LFs, shown in Fig. [l6| are found to have different shape s than the 
central LFs (solid lines). Mobasher & Trentham (19981 have de- 



rived a K band LF for the field around NGC 4839. However, due 
to a small field size (9.9 arcmin 2 ) and large uncertainties in back- 
gro und su btraction, their LF is essentially unconstrained. Lobo et 
al. (1997) have derived a V band LF with a faint end slope compa- 
rable to their central LF, using an area of 117 arcmin 2 . We find a 
B band LF which is much shallower than the central LF. The faint 
end of the U band LF is also flatter, except for the last point. In the 
r band we have a lack of galaxies fainter than Mr = — 16, but the 
faint end slope is comparable to the slope of the central LF. 



9 DISCUSSION 

In this paper we have used a statistical method to determine in U, B 
and r the LF of the galaxies in the Coma cluster and the dependence 
of the LF on projected distance from the cluster centre. We did find 
changes in shape as a function of radius indicating changes in the 
galaxy population in the cluster. At the same time this implies that 
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differences. 



5 log h 100 M B - 5 log h 100 U u - 5 log h 100 

for the area around NGC 4839. The best-fitting Schechter functions of the central luminosity functions are overplotted to emphasize shape 



comparison with other work is difficult as the LF shape depends 
critically on the area chosen for study. The changes in the LF shapes 
are very likely related to the local environment. We will accordingly 
discuss this possibility in more detail here. 

The a parameter of the Schechter function measures effec- 
tively the faint end slope of the luminosity distribution and charac- 
terizes the intermediate and dwarf galaxy populations. A flat faint 
end slope i mplies a lack of these low luminosity galaxies. Lopez- 
Cruz et al. ( 1997) observe the trend that steep faint end slopes are 



detected in poorer clusters and the flatter slopes are, on average, 
found in richer clusters. This suggests that environmental proper- 
ties could dictate the faint galaxy population. In a scenario in which 
mergers, interactions, tidal stripping, destruction of dwarfs, infall 
etc. play a role we expect to see a lack of faint galaxies towards the 
regions of higher galaxy density. This then should be reflected in a 
flattening of the low luminosity end of the LF. We do observe this 
effect; e.g. the slope of the faint end of the U band LF decreases 
from a = 0.20 at ~ 0.7 /i^oo Mpc to a ~ 0-13 in the centre. 
In terms of the slope of the Schechter function these values corre- 
spond to a = —1.5 and a = —1.32. The effect is also present 
in the B and r band LFs, and measured out to a larger radius of 
r = 1.3 fc^oo Mpc. 

When comparing the faint end slopes of the Coma LFs with 
those from the SDSS the most striking result is that the r and B 
band LFs are very similar. This is remarkable given the differ- 
ences in environment and galaxy populations. However, the U band 
slopes of Coma are steeper than in the field so apparently there is 
a relation between the relative increase in low luminosity systems 
in Coma and the colour of the band. A possible interpretation is 
that the (infalling?) dwarf galaxies in the outer parts of Coma are 
undergoing bursts of star formation triggered by interactions with 
neighbour galaxies and/or the intra-cluster medium. As a result the 
galaxies brighten, preferentially in the U band, causing a steepen- 
ing of the faint end slopes of the LFs. If correct, this is a clear sign 
that the Coma cluster is not a relaxed system, but that, especially in 
its periphery, the cluster is still forming and inducing strong evo- 



lution to the galaxy population. The flattening of all Coma LFs to- 
wards the cluster centre hints that there the galaxies have already 
lost most of their gas and enhanced star formation has long ceased. 

The group around NGC 4839 has been subject of debate in 
the literature, because it is not clear whether it is falling into the 
cluster for the first time or has already made one pass through (for 
scenarios see e.g. Colless & D unn 19 96; White et al. 1993; Bums et 
al. 1994). Bravo- Alfaro et al. (200C) have done Hi imaging of the 
Coma cluster and the NGC 4839 group. From the nondetections in 
the close vicinity of NGC 4839 and the presence of several starburts 
and post-starbust galaxies with very low Hi content in that zone 
they conclude that it passed at least once through the core. The LF 
of the field around NGC 4839 has a different faint end slope than 
the central LFs which could be related to the dynamical history of 
this group. The observed shape differences could be explained if 
during the passage a large fraction of the dwarf galaxies has been 
stripped from the group and redistributed throughout the cluster 
potential. 

To explain the dependence of the shape of the LF on pro- 
jected distance from the cluster centre we would need much more 
information. Knowledge of the galaxy morphologies, redshifts, 
the morphology-density relation (Dressier 1980; Whitmore et al. 
1983), the type-dependent LFs, Hi observations etc. are necessary 
to derive a complete scenario in which mergers, tidal stripping, in- 
fall etc. play a role. For instance, the observed dips and increasing 
faint end slopes could be the combined result of Coma's morpho- 
logical composition together with the shapes of the type-dependent 
LFs. 



10 SUMMARY AND CONCLUSIONS 

We have presented the first results of a wide field photometric sur- 
vey of the Coma cluster in the U, B and r bands. The derivation of 
the source catalogue, along with the steps concerning the pipeline 
reduction, have been discussed. Our high quality data provides a 
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valuable low z comparison sample for studies of galaxy morphol- 
ogy, colour and luminosity at higher z. 

In this paper we have used the data to study the dependence of 
the galaxy LF on passband and projected distance from the cluster 
centre. The LFs of the complete data set cannot be represented by 
single Schechter functions. The central U, B and r band LFs can 
be represented by Schechter functions with parameters as listed in 
Table ^ The expectation that the shape of the LF depends on envi- 
ronment is confirmed. The LFs as a function of distance from the 
cluster centre have (very) different shapes than the central LFs and, 
therefore, no universal general LF exists. The faint ends of the LFs 
become steeper towards the outskirts of the cluster. A steepening of 
the faint ends of the LFs towards less dense regions clearly supports 
the existence of environmental effects. The difference in faint end 
slopes of the overall LFs and those of the field is colour dependent 
(strongest in U) and can be attributed to enhanced star formation in 
the dwarf galaxy population in the outer parts of Coma where evo- 
lution apparently is still very strong. The LFs of the field around 
NGC 4839 support the idea that this group has passed through the 
cluster centre at least once. 

We are in the process of creating a catalogue of all spectro- 
scopically confirmed cluster members. This will enable us to revise 
the Coma LFs and to check the quality of the statistical method to 
determine LFs. We have also obtained Westerbork Synthesis Ra- 
dio Telescope (WSRT) Hi mosaic data covering an total area of 
2.94° x 2.05° in order to study the Hi properties as function of 
environment and assess the importance of merging and stripping. 
Combined with our optical data set this will greatly enhance the 
ability to study the structure and dynamics of Coma, using a data 
set that is unrivalled by what is available for any other cluster. 



ACKNOWLEDGEMENTS 

We thank the referee for detailed and constructive comments. Part 
of the data was made publically available through the Isaac New- 
ton Groups' Wide Field Camera Survey Programme. P. G. v. D. 
acknowledges support by NASA through Hubble Fellowship grant 
HF-01 126. 01-99 A awarded by the Space Telescope Science Insti- 
tute, which is operated by the Association of Universities for Re- 
search in Astronomy, Inc., for NASA under contract NAS 5-26555. 



Dressier A., 1980, ApJ, 236, 351 

Driver S.P., Couch W.J., Phillipps S., 1998, MNRAS, 301, 369 

Fukugita M., Ichikawa T., Gunn J.E., Doi M., Shimasaku K., Schneider 

D. P., 1996, AJ, 111, 1748 
Geller M.J., et al., 1997, AJ, 1 14, 2205 

Godwin J. G., Metcalfe N., Peach J. V., 1983, MNRAS, 202, 1 13 
Landolt A.U., 1992, AJ, 104, 340 

Lobo C, Biviano A., Durret R, Gerbal D., Le Fevre O., Mazure A., Slezak 

E. , 1997, A&A, 317, 385 

Lopez-Cruz O., Yee H.K.C., Brown J.P., Jones C. Forman W., 1997, ApJ, 

475, L97 
Lin H. et al., 1996, ApJ, 464, 60 

Lin H., Yee H.K.C., Carlberg R.G., Ellingson E., 1997, ApJ, 475, 494 
Lin H., Yee H.K.C., Carlberg R.G., Morris S.L., Sawicki M., Patton D.R., 

Wirth G., Shepherd C.W., 1999, ApJ, 518, 533 
Loveday J., Peterson B.A., Efstathiou G, Maddox S.J., 1992, ApJ, 390, 338 
Lugger P.M., 1989, ApJ, 343, 572 

Marzke R.O., Nicolaci da Costa L., Pellegrini PS., Willmer C.N.A., Geller 
M.J., 1998, ApJ, 503, 617 

Mazure A., Casoli F, Durret F, Gerbal D., eds., 1998, "Untangling Coma 
Berenices: A New Vision of an Old Cluster", proceedings of the meet- 
ing held in Marseille (France), June 17-20, 1997, World Scientific Pub- 
lishing Co. Pte. Ltd. 

McGaugh S.S., Bothun G.D., 1994, AJ, 107, 530 

Mobasher B., Trentham N.D., 1998, MNRAS, 293, 315 

Pickering T.E., Impey CD., van Gorkom J.H., Bothun G.D., 1997, AJ, 1 14, 
1858 

Schechter P., 1976, ApJ, 203, 297 

Schlegel D.J., Finkbeiner D.P., Davis M., 1998, ApJ, 500, 525 
Seeker J., Harris W.E., Plummer J.D., 1997, PASP, 109, 1377 
Sprayberry D., Impey CD., Bothun G.D., Irwin M.J., 1995, AJ, 109, 558 
Thompson L.A., Gregory S.A., 1993, AJ, 106, 2197 
Trentham N.D., 1998a, MNRAS, 293, 71 
Trentham N.D., 1998b, MNRAS, 294, 193 

van Dokkum P.G., Franx M, Kelson D.D., Illingworth G.D., Fisher D., Fab- 

ricant D., 1998, ApJ, 500, 714 
White S.D.M., Briel U.G., Henry J.P., 1993, MNRAS, 261, L8 
Whitmore B.C., Gilmore D.M., Jones C, 1993, ApJ, 407, 489 
Wolf M., 1901, Astron. Nachr., 155, 127 



REFERENCES 

Andreon S., 1999, A&A, 351, 65 
Andreon S., Pello R., 2000, A&A, 353, 479 
Andreon S., 1998, A&A, 336, 98 

Bernstein CM., Nichol R.C., Tyson J.A., Ulmer M.P., Wittman D., 1995, 
AJ, 110, 1507 

Biviano A., Durret F, Gerbal D., le Fevre O., Lobo C, Mazure A., Slezak 

E., 1995, A&A, 297, 610 
Bravo-Alfaro H., Cayatte V., van Gorkom J.H., Balkowski C, 2000, AJ, 

119, 580 

Bertin E., Arnouts S., 1996, A&AS, 117, 393 

Blanton M.R., et al., 2001, AJ, 121, 2358 

Bower R.C, Lucey J.R., Ellis R.S., 1992, MNRAS, 254, 589 

Burns J.O., Roettiger K., Ledlow M., Klypin A., 1994, ApJ, 427, L87 

Colless M., Dunn A. M., 1996, ApJ, 458, 435 

de Blok W.J.G., van der Hulst J.M., Bothun CD., 1995, MNRAS, 274, 235 
de Propris R., Eisenhardt PR., Stanford S.A., Dickinson M., 1998, AJ, 503, 
L45 

Deul E. R., 1998, Pipeline Documentation Ver. 1.1, Leiden Observatory 
de Vaucouleurs G, 1948, Ann. d'Astrophys., 11, 247 



© 2001 RAS, MNRAS 000, 000-000 



